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We study the (K~,p) reaction on 12 C with a kaon beam of 1 GeV momentum, paying 
a special attention to the region of emitted protons having kinetic energy above 600 MeV, 
which was used to claim a deep kaon nucleus optical potential [l[ • The experiment looks for 
fast protons emitted from the absorption of in flight kaons by nuclei, but in coincidence with 
at least one chargedparticle in the decay counters sandwiching the target. The analysis of 
the data is done in [l| assuming that the coincidence requirement does not change the shape 
of the final spectra. However our detailed calculations show that this assumption doesn't 
hold, and, thus, the final conclusion of this experiment is doubtful. 

We perform Monte Carlo simulation of this reaction. The advantage of our method with 
respect to Green's function method used in [l[ is that it allows to account not only for 
quasi-elastic K~p scattering, but also for the other processes which contribute to the proton 
spectra. We investigated the effect of the multi-scatterings and of the K~ absorptions by one 
and two nucleons (K~N — > irY and K~NN — * YN) followed by the decay of the hyperon 
in irN. We show that all these mechanisms allow us to explain reasonably well the observed 
spectrum with standard shallow kaon nucleus optical potential, obtained in chiral models. 

PACS numbers: 13.75.-n,12.39.Fe,14.20.Jn,11.30.Hv 
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The issue of the kaon interaction in the nu- 
cleus has attracted much attention in past years. 
Although from the study of kaon atoms one 
knows that the -fT - -nucleus potential is attrac- 
tive [2j, the discussion centers on how attractive 
the potential is and whether it can accommodate 
deeply bound kaon atoms (kaonic nuclei), which 
could be observed in direct reactions. 

All modern potentials based on underlying 
chiral dynamics of the K N interaction 0, 3, la, 
H, 0] lead to moderate potentials of the order 
of 60 MeV attraction at normal nuclear density. 
They also have a large imaginary part making 
the width of the bound states much larger than 
the energy separation between the levels, which 
would rule out the experimental observation of 
these states. 

Deep K~-N optical potentials are preferred 



by the phenomenological fits to kaon atoms data. 
One of the most known extreme cases of these 
type is a highly attractive phenomenological po- 
tential with about 600 MeV strength in the center 
of the nucleus, introduced in [8|, |9j. In these pic- 
ture such an attractive K~, inserted inside the 
nucleus, would lead to a shrinkage of the nucleus, 
generating a new very compact object - kaonic nu- 
cleus - with central density which can be 10 times 
larger than normal nuclear density. Such super- 
deep potentials were criticized in Tol. [Til. Il2l. flU] . 

It is important to keep in mind that in kaon 
atoms the K~ is primary bound by the Coulomb 
force. These are extended systems, and therefore 
their properties cannot directly tell us about the 
-ftT~-N potential at short distances. From the ex- 
perimental side the search for bound K~ states 
with nucleons is a most direct and clear way to 
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answer whether the isT~-nucleon potential is deep 
or shallow, because only deep potential may gen- 
erate states sufficiently narrow to be observed ex- 
perimentally. Experimental attempts to resolve 
this situation are made since 2004, but the situ- 
ation is still very unclear. 

Several claims of observed deeply bound K~ 
states have been made. However the first one, 
about K~pnn state bound by 195 MeV from the 
experiment at KEK 14| , is now withdrawn after a 
new more precise experiment [15]. The peaks seen 
by FINUDA and originally interpreted in terms of 
deeply bound K~pp [l(| and K~ppn [TtJ clusters, 
are now pu t under the question, because in Refs. 
HI, |2(J Hi] these peaks found explanations 



H 

based on conventional reactions that unavoidably 
occur in the process of kaon absorption. 

There are also claims (with very low sta- 
tistical significance) of K~pp and K~ppn bound 
states from p annihilation in 4 He at rest mea- 
sured by OBELIX@CERN Q. And the recent 
claim of K~pp bound state, seen in pp — > K + X 
reaction, from DISTO experiment [23|]. These ex- 
perimental claims are under investigation now. 
Before calling in new physics one has to make 
sure that these data cannot be explained with 
conventional mechanisms. 

There is however one more experiment 
where the authors claim the evidence for a strong 
kaon-nucleons potential, with a depth of the or- 
der of 200 MeV [1!]. The experiment looks for 
fast protons emitted from the absorption of in 
flight kaons by 12 C in coincidence with at least 
one charge particle in the decay counters sand- 
wiching the target. The data analysis in [l(] is 
based on the assumption that the coincidence re- 
quirement does not change the shape of the fi- 
nal spectra. We shall see that this assumption 
doesn't hold and the interpretation of the data 
requires a more thorough approach than the one 
used in that work. 

One of the shortcomings of Ref. [l[ stems 



from employing the Green's function method 24] 
to analyze the data. The only mechanism con- 
sidered in Ref. [l[ for the emission of fast protons 
is the Kp — ► Kp process, taking into account the 
optical potential for the slow kaon in the final 
state. We shall show that there are other mech- 
anisms that contribute to generate fast protons, 
namely multi-scattering reactions, and kaon ab- 
sorption by one nucleon, K~N — > 7rS or K~N — > 



7rA or by a pair of nucleons, KNN — ► EiV and 
KNN — > AN, followed by decay of E or A into 
ttN. The contributions from these processes were 
also suggested in Ref. [HI]. 

In the present work, we take into account 
all the above mentioned reactions by means of 
a Monte Carlo simulation [2f|. The election of 
which reaction occurs at a certain point in the 
nucleus is done as usual. One chooses a step size 
51 and calculates, by means of Oip5l, the proba- 
bilities that any of the possible reactions happens 
i = Quasi- elastic, IN absorption, 2N absorption; 
p is nucleon density. The size of 51 is small enough 
such that the sum of probabilities that any reac- 
tion occurs is reasonably smaller than unity. A 
random number from to 1 is generated and a 
reaction occurs if the number falls within the cor- 
responding segment of length given by its proba- 
bility. If the random number falls outside the sum 
of all segments then this means that no reaction 
has taken place and the kaon is allowed to pro- 
ceed one further step 51. The simulation of one 
event is over when all the produced particles leave 
the nucleus. To adapt the calculations to the ex- 
periment of [l| we select "good events" with fast 
protons that emerge within an angle of 4.1 de- 
grees in the nuclear rest frame (lab frame). As in 
Pi we plot our obtained l2 C(K~ ,p) spectrum as 
a function of a binding energy of the kaon, Eg, 
should the process correspond to the trapping of 
a kaon in a bound state and emission of the fast 
proton. 

If there is a quasi-elastic collision at a certain 
point, then the momentum of the K~ and that of 
the nucleon, which is randomly chosen within the 
Fermi boosted to their CM frame. The 

direction of the scattered momenta is determined 
according to the experimental cross section. A 
boost to the lab frame determines the final kaon 
and nucleon momenta. The event is kept as long 
as the size of the nucleon momentum is larger 
than the local Fermi momentum. Since we take 
into account secondary collisions we also consider 
the reaction K~p — > K°n and K~n — » K~n with 
their corresponding cross sections. 

Once primary nucleons are produced they 
are also followed through the nucleus taking into 
account the probability that they collide with 
other nucleons, losing en ergy and changing their 
direction, see [181 . 1 19 . |2CJ, |2l| for more details. 

We also follow the rescattered kaon on its 
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way through the nucleus. In the subsequent 
interaction process we let the kaon experience 
whichever reaction of the three that we con- 
sider (quasi-elastic, one-body absorption, two- 
body absorption) according to their probabilities. 
This procedure continues until kaon is absorbed 
or leaves of the nucleus. 

Apart from following the kaons and nucle- 
ons, our calculations also need to consider the 
quasi-elastic scattering of A's and E's (produced 
in the kaon absorption reactions) on their way 
through the residual nucleus. Given the uncer- 
tainties in the hyperon-nucleon cross sections, we 
may use for XiV scattering the relation cr^N = 
2cjvjv/3, based on a simple non-strange quark 
counting rule. In the case of AiV scattering, we 
use the refined parameterization of Ref. [27], as 
was also done in Ref. [20]. 

One nucleon K~ absorption leads to K~N 
— * 7rA or K~N — > 7rS, with all the possible 
charge combinations. The elastic and inelastic 
two-body KN cross sections for kaons are taken 
from the Particle Data Group [28| . 

The kaon absorption by two nucleons is a bit 
more tricky. Here we take into account the follow- 
ing processes: K^NN -> AiV or K'NN -» T,N 
with all possible charge combinations. In these 
reactions an energetic nucleon is produced, as 
well as a A or a S. Both the nucleon and the 
hyperon are followed through the nucleus as dis- 
cussed above. Once out of the nucleus, the hyper- 
ons are let to decay weakly into ttN pairs. There- 
fore, the two-body absorption process provides 
a double source of fast protons, those directly 
produced in two nucleon absorption reaction and 
those coming from hyperon decays. 

We assume a total two body absorption rate 
to be 20% that of one body absorption at about 
nuclear matter density, something that one can 
infer from data of K~ absorption in 4 He [29] . In 
practice, this is implemented in the following way. 
The probability per unit length for two nucleon 
absorption is proportional to the square of the 
nucleon density: Pk-nn(p) = Cabs/0 2 . We as- 
sume that (hk-nn) = C ab s(p 2 ) = 0.2(p K - N ) = 
0.2cr^l N (p) , where c^ljy accounts for the to- 
tal one nucleon absorption cross section and, in 
symmetric nuclear matter it is given by: c^l^ = 

^ at K-p Jr ° X K-n~ a K-p^K-p — a K-n^K~n)/2 ■ Tak- 
ing (p) = po/2, where po = 0.17 fm~ 3 is normal 



nuclear matter density, we obtain C a b s ~ 6 fm . 

The different partial processes that can 
take place in a two-nucleon absorption reac- 
tion are: K~pp — > pA, K~pn — ► 
nA, pYi~; K~nn — > n£~ . Ideally, their cor- 
responding branching ratios should be obtained 
from relevant microscopic mechanisms, however 
in the present exploratory work, we will consider 
a much simpler approach consisting of assigning 
equal probability to each of the above reactions. 
Noting that the chance of the kaon to find a pn 
pair is twice as large as that for pp or nn pairs, 
we finally assign a probability of 3/10 for having 
a pT, pair in the final state, 4/10 for nS, 1/10 for 
pA and 2/10 for nA. 
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FIG. 1: Calculated 12 C{K~,p) spectra with V opt = 
(—60, — 60)p/po MeV, taking into account only quasi- 
elastic processes (dash-dotted line), and including all 
the contributing processes (full line). 

We also take into account a kaon optical po- 
tential V op t = Re Vopt + i Im V pt , which will influ- 
ence the kaon propagation through the nucleus, 
especially when it will acquire a relatively low mo- 
mentum after a high momentum transfer quasi- 
elastic collision. In the present study we take the 
strength of the potential as predicted by chiral 
models: ReV opt = -60p/p MeV 0, II & 0, H 5 
ImVopt ~ — 60p/po MeV, as in the experimental 
paper [l| and the theoretical study of jj|. 

In the Monte Carlo simulation we implement 
this distribution by generating a random kaon 
mass Mk around a central value, Mk + Re Vopt, 
within a certain extension determined by the 
width of the distribution Tk = — 2ImV pt- The 
probability assigned to each value of Mk follows 
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the Breit-Wigner distribution given by the kaon 
spectral function: 
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In Fig. [T] we show the results of the Monte 
Carlo simulation obtained with an optical po- 
tential Kpt = (— 60, — 60)p/po MeV: first, tak- 
ing into account only quasi-elastic processes; and 
then taking into account all the discussed mech- 
anisms. We can see that there is some strength 
gained in the region of " bound kaons" due to the 
new mechanisms. Although not shown separately 
in the figure, we have observed that one nucleon 
absorption and several rescatterings contribute to 
the region — Eb > —50 MeV. To some extent, 
this strength can be simulated by the parametric 
background used in However, this is not true 
anymore for the two nucleon absorption process, 
which contributes to all values of —Eb, starting 
from almost as low as —300 MeV. 

It is very important to keep in mind that 
in the spectrum of [l(] the outgoing forward pro- 
tons were measured in coincidence with at least 
one charged particle in the decay counters sand- 
wiching the target. Obviously, the real simula- 
tion of such a coincidence experiment is tremen- 
dously difficult, practically impossible with high 
accuracy, because it would require tracing out all 
the charged particles coming out from all possible 
scatterings and decays. Although we are study- 
ing many processes and following many particles 
in our Monte Carlo simulation, which is not the 
case in the Green function method used in the 
data analysis [If , we cannot simulate precisely the 
real coincidence effect. 

The best we can do is to eliminate the pro- 
cesses which, for sure, will not produce a coinci- 
dence, this can be called minimal coincidence re- 
quirement. If the kaon in the first quasi-elastic 
scattering produces an energetic proton falling 
into the peaked region of the spectra, then the 
emerging kaon will be scattered backwards. In 
our Monte Carlo simulations we can select events 
were neither the proton, nor the kaon will have 
any further reaction after such a scattering. In 
these cases, although there is a "good" outgoing 
proton, there are no charged particles going out 
with the right angle with respect to the beam axis 
to hit a decay counter, since the K~ escapes un- 
detected in the backward direction. Therefore, 
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FIG. 2: Calculated 12 C{K~,p) spectra with V opt = 
(—60, — 60)p/po MeV taking into account all con- 
tributing processes (dash-dotted line). Then we 
impose minimal coincidence requirement (full line). 
Data points are from [l| . 



this type of events must be eliminated for com- 
parison with the experimental spectra. 

It is clear from Fig. Q] that the main source 
of the energetic protons for 12 C(K~ ,p) spectra 
is K~p quasi-elastic scattering, however many of 
these events will not pass the coincidence condi- 
tion. Implementing the minimal coincidence re- 
quirement, as discussed above, we will cut off a 
substantial part of the potentially "good" events, 
and drastically change the form of the final spec- 
trum, as illustrated in Fig. [2j 

To further simulate the coincidence require- 
ment we introduce additional constant suppres- 
sion factors to the obtained spectrum - see Fig. 
[3l Comparing our results with the experimen- 
tal data we can conclude that in the "bound" 
region, — Eb < MeV, these additional suppres- 
sion is about ~ 0.7 and more or less homogeneous, 
while in the continuum the suppression weakens 
and for — Eb > 50 MeV it is negligible. This pic- 
ture is natural from the physical point of view, 
because the r.h.s. of the spectrum, Fig. [31 with 
relatively low momentum protons is mostly pop- 
ulated by many particle final states, which have 
a good chance to score the coincidence. 

To conclude, the main point of our analysis 
is not to state that the data of Ref. [l[ supports 
Re Vopt = -60p/Vo MeV rather than -200p/p . 
We want to make it clear that trying to simulate 
these data one necessarily introduces large un- 
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FIG. 3: Calculated 12 C(K~,p) spectrum with V opt = 
(—60, — 60)p/po MeV with minimal coincidence re- 
quirement - solid line; and with additional suppression 
factors - dash-dotted and dotted lines. Experimental 
points are from [l[ ■ 



certainties due to the experimental set up. Thus, 
this experiment is not appropriate for extracting 
information on the kaon optical potential. 

Contrary to what it is assumed in Ref. [l|, 
we clearly see, Fig. O that the spectrum shape is 
affected by the required coincidence. The experi- 
mental data without the coincidence requirement 
would be a more useful observable. 
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